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On the basis of Boltzmann's equation, and including anisotropic scattering in the collision opera- 
tor, we investigate the effect of one-dimensional superlattices on two-dimensional electron systems. 
In addition to superlattices defined by static electric and magnetic fields, we consider mobility su- 
perlattices describing a spatially modulated density of scattering centers. We prove that magnetic 
and electric superlattices in ^-direction affect only the resistivity component p xx if the mobility 
is homogeneous, whereas a mobility lattice in x-direction in the absence of electric and magnetic 
modulations affects only p yy . Solving Boltzmann's equation numerically, we calculate the positive 
magnetoresistance in weak magnetic fields and the Weiss oscillations in stronger fields within a 
unified approach. 



I. INTRODUCTION 

Two-dimensional electron systems (2DESs) in high- 
mobility Ali_ x Ga x As-GaAs heterostructures with one- 
dimensional lateral superlattices show interesting magne- 
totransport properties at low temperatures (T « 4K or 
less) . Experiments on typical samples with electron 
density n e i ~ 3 • 10 11 cm~ 2 , mobility ji ~ 10 6 cm 2 /Vs, 
and a weak electric modulation JlJ in a;-direction of pe- 
riod a ~ 300 nm revealed the following features. If the 
homogeneous magnetic field Bq applied perpendicular to 
the 2DES is very weak (Bq < 0.03T), one finds a pro- 
nounced positive magnetoresistance in the component 
Pxx- With increasing Bq, this saturates and is followed 
by the "Weiss oscillations" (0.1T < B < IT), which at 
higher fields (typically B > 0.6T at T = 4.2K) are su- 
perimposed by the Shubnikov-de Haas oscillations. The 
resistivity component p yy does not exhibit such a pos- 
itive magnetoresistance, but shows Weiss oscillations of 
opposite phase. No significant effect of the modulation 
on the Hall resistance has been observed. Experiments 
on samples with a stronger electric modulation ||,Q , and 
also the recent experiments on magnetically modulated 
systems H|), led to similar results, although they usually 
did not provide the full information about the conductiv- 
ity tensor. The Weiss oscillations were first understood 
within a quantum mechanical picture : The modula- 
tion lifts the degeneracy of the Landau levels and leads to 
dispersive Landau bands of oscillatory widths. The group 
velocity of these bands leads to a "band conductivity" , 
which explains the Weiss oscillations of p xx . The oscilla- 
tory density of states affects the scattering rate and gives 
rise to the Weiss oscillations of p vy (?|!| . Beenakker |l(]] 
pointed out that the Weiss oscillations of p xx can also be 
understood classically in terms of a guiding center drift 
of cyclotron orbits in the electric modulation field. He 
proved from Boltzmann's equation that, under the as- 
sumption of a constant relaxation time, p xx is the only 
component of the resistivity tensor that is affected by 
the modulation. The low-i?o positive magnetoresistance 



has also been explained classically 0, as being due to 
channeled trajectories which, in x-direction are localized 
within a single period of the modulation. Beenakker ne- 
glected these channeled orbits in his explicit calculation 
and could therefore not explain the positive magnetore- 
sistance, whereas Beton et al. Q could not extend their 
calculation to the regime of Weiss oscillations. So far 
only the semi-classical approach by Stfeda and cowork- 
ers — p~3|] has been applied to both regimes. However, 
this approach relies on quantum concepts such as en- 
ergy band structure and magnetic breakdown. A unified 
calculation, which explains both the low-Bo magnetore- 
sistance and the Weiss oscillations of p xx totally within 
classical dynamics, is still missing. The purpose of this 
paper is to present such a classical calculation. 

We will extend Beenakker's work by including 
other modulation sources in addition to a periodic elec- 
tric field. We will consider a magnetic modulation and 
also a "mobility modulation" , i.e. a periodic position 
dependence of the collision operator C of Boltzmann's 
equation, describing a spatial modulation of the density 
of scatterers. We will also allow for anisotropic scatter- 
ing of the electrons by, e.g., impurity potentials of finite 
range, and thus go beyond the relaxation time approx- 
imation of Boltzmann's equation. Within our classical 
approach we will prove that, in the absence of a mobility 
modulation, a one-dimensional superlattice defined by x- 
dependent electric and magnetic fields affects only p xx , 
even if anisotropic scattering is included. Allowing for a 
periodic mobility modulation in x-direction, we can ob- 
tain oscillations of p yy , however with the same phase as 
those of p xx in the case of a pure electric modulation. 

After giving some details about our model in Sect. II, 
we present the formal calculations revealing the struc- 
ture of the magnetoresistance tensor in Sect. Ill and 
in Appendix [a|. In Sect. IV we present analytical re- 
sults: For the regime of Weiss oscillations, we generalize 
Beenakker's calculation to include magnetic and mobility 
modulations, and compare with simplified calculations 
based on the evaluation of the drift velocities of cyclotron 
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orbits @|3jT5). For the regime of channeled trajecto- 
ries we concentrate on the collisionless limit, where we 
obtain an algebraic dependence of the conductance on 
the modulation strength. Therefore, in contrast to the 
regime of Weiss oscillations, a simple power expansion of 
the transport coefficients with respect to the modulation 
strength is not possible in this regime. 

Numerical results based on a Fourier expansion of 
Boltzmann's equation are presented in Sect. V. For suf- 
ficiently weak modulation strengths, these results cover 
the whole range of low-field magnetoresistance and Weiss 
oscillations, and they include the effect of anisotropic 
scattering as well as of electric, magnetic and mobility 
modulations. In a final Sect. VI we summarize the most 
important results. 



II. THE MODEL 

It is known from quantum mechanical calculations 
|j],D that the Weiss oscillations on typical high-mobility 
samples are observed in a peculiar temperature range. 
The temperature must be high enough, so that in the 
regime of Weiss oscillations the Shubnikov-de Hass os- 
cillations are not resolved. On the other hand, it must 
be low enough, so that the transport coefficients are de- 
termined only by the electronic states near the Fermi 
energy [p| Jl^Jl5t| . To describe the relevant physics in this 
temperature window, we consider the 2DES in the x- 
y plane as a degenerate Fermi gas, with Fermi energy 
E F = (m/2)v 2 F , of ( non-interacting) particles with effec- 
tive mass m and charge — e obeying classical dynamics. 
The velocity v — r = (x,y, 0) of an electron obeys New- 
ton's equation 



mv = -e [F + (v x B) /c] . 



(2.1) 



To describe a modulated system with period a in equi- 
librium, we write the electric field as F(x) = W(x)/e 
and the magnetic field as B(x) = (0,0, Bo + B m (x)), 
and assume that the periodic functions V{x) = V(x + a) 
and B m (x) = B m (x + a) have zero average values. The 
thermal equilibrium at constant temperature T and elec- 
trochemical potential fi* is described by the distribution 
function 



/ eg (r,v)=/ (£;(r ) v);T, / i*), 



(2.2) 



where the Fermi function fo(E; T, ^*) = 1/[1 + exp[(E — 
H*)/kBT}] depends on the dynamical variables only 
through the conserved energy 



E(r,x) = (m/2)v 2 + V(x) . 



(2.3) 



In more general situations, the distribution function 
/(r, v, t), which yields electron density n e i and current 
density j according to 

Ttej(r,t) = (/>« , j(r,i) = (-ev/>„ , (2.4) 



where (f) v = \2m 2 j [2-nK) 2 \ J d 2 vf, is determined by 
Boltzmann's equation 



and suitable boundary conditions. For any reasonable 
microscopic model of elastic and inelastic scattering [[l6| , 
the collision operator C vanishes after averaging over 
the velocity, (C% = 0. As a consequence, the v- 
integral over Eq. (2.5) yields the continuity equation, 
—edn e i/dt+\/ ■ j = 0. If we want to choose a phenomeno- 
logical model for C, we should be careful not to destroy 
this property. 

In the following we consider the stationary, linear re- 
sponse of the modulated system to an external homoge- 
neous electric field E° , which describes the average effect 
of a voltage applied across the sample. The result will be 
a non-zero current density j(x) and a correction Sn e i(x) 
to the equilibrium electron density. At low temperatures 
[ksT -Cu* — V"( x)1, the equilibrium density follows from 
Eqs. (2Jj) and (2.4) as n e i(x) = D [E F - V(x)], where 
Dq = m/ (nh 2 ) is the density of states and Ep = fj,*(T — > 
0) the Fermi energy of the homogeneous 2DES. 

With the ansatz 



Wr, v) = f Q (E(r, v); T, //) - e00(r, v) (2.6) 

for the stationary solution, and the assumption that the 
correction 0(r, v) to the equilibrium distribution is linear 
in E , we obtain the linearized Boltzmann equation 



%v)-C[^;r,v]=vE° ; 



(2.7) 



where T> 



■ d/dr + v • d/dv, and v is determined by 
Eq. (2.1) with the force in the equilibrium state. For 
T — > 0, dfo/d^i* -> 5(E(r,v) - E F ) and with Eq. (|J), 
which implies for the magnitude of the vel ocity v{x) — 
v F [l-V(x)/E F } 1 / 2 , the 2D w-integral in Eq. (§!]) reduces 
to an integral over the polar angle in the velocity space, 
e.g. 



e 2 D 



^■v(x)vL(<p)<f>(x,<p), (2.8) 



where u((p) — (cos ip, sin if), and the dependence of 4> on 



Ep is not explicitly indicated. The drift term of Eq. (2.7) 
reads in these variables 

9 9 
V = v(x) cos ip— + [oj c + uj m (x) +ui e (x,(f)] — , (2.9) 
ox dip 

with co c = eBo/mc the cyclotron frequency due to the 
average magnetic field, and ui m (x) — eB m (x)/mc and 
u> e (x, tp) = — shxipdv/dx are due to the magnetic and 
the electric modulation, respectively. The collision oper- 
ator describes the change of the distribution function due 
to scattering events with a spatial extent which is short 
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on the scale of the drift motion Jig] . If there are dif- 
ferent scattering mechanisms, we may write C[<p;x,p] = 
J2jC 3 [(j);x,(p}, with 

Cj[<f>;x,ip\ = -^j-^ I — Pj(ip' - tp)[<t>(x,tp') - (f)(x,ip)] 



]Tt| , P~8|] . Neglecting <f>i oc would describe relaxatio n t o- 



Tj(x) J 2tt 



(2.10) 



where we have expressed the renormalized differential 
cross sections in terms of relaxation times Tj (x) and di- 
mensionless kernels Pj(p)- Here, the last term describes 
the decay of the distribution function due to scattering 
processes of the j-th type which start in phase space 
at (x, ip), whereas 4>(x' , ip') describes the back-scattering 
into this state. Assuming that the distribution function 
changes little on the range of the scattering processes, 
we consider only x' = x in the back-scattering term. For 
scattering by impurities, the scattering rate is propor- 
tional to the impurity density, which may be a function 
of the position coordinate x in periodically modulated 
samples. In order to describe a periodic modulation of 
the electron mobility, we write 1/Tj(x) — l/fj + Tj{x) 
and assume that the oscillating parts rj (x) — rj (x + a) 
of the scattering rates have zero average values. Further 
we assume that the kernels Pj(p), with the normalization 
S-TrdfPj{(p) = 2tt, have the symmetry Pj(— (p) = Pj((p). 
This follows from the microscopic symmetry if the effect 
of the magnetic field during an individual scattering pro- 
cess is neglected. 

It is interesting to note that the trigonometric func- 
tions are eigenfunctions of the collision operator: for 
VVnCv 9 ) = exp(imip) one has 



C[ip m ;x,(p] = -[l/r^\x)} ip m (tp) , 
with the eigenvalues 



iAi m) (*) 



where = (2tt)^ 1 dtpPj(tp) coa(mp). The same 
holds for the real part and the imaginary part of tp m 
separately. We will demonstrate in the following that 
several transport properties do not depend on all th e de- 
tails of the collision operator defined by Eq. ( 2.10 ), but 



(2.11) 



(2.12) 



only on the transport time T tr (x) = t^\x), which we 
write in the form 1/T tr (x) = l/f tr + r tr (ir). For instance, 
the Drude conductivity tensor of the homogeneous 2DES, 
with rtr(x) = 0, depends only on ft r - Thus, the trans- 
port scattering rates due to the different scattering mech- 
anisms simply add. 

The special case Pj(p) = 1 describes isotropic scat- 
tering, and is identical with the relaxation time approx- 



imation, Cj [4>] = 



4>io C )/Tj, where the relaxation 



is towards the local equilibrium distribution <pi oc (x) 
J" d<p<f)(x, <p)/2tt. This is the correct form of the re- 
laxation time approximation for inhomogeneous systems 



wards the total equilibrium distribution [cf. Eq. (2.6)], 
and would violate explicitly the equation of continuity. 
As we will discuss below, this can lead to qualitatively 
wrong results for the tranport coefficients. 



III. STRUCTURE OF THE RESISTIVITY 
TENSOR 

In this section we prove from Boltzmann's equation 
that a combined electric and magnetic modulation in in- 
direction affects only the component p xx of the effective 
resistivity tensor, if there is no position dependence of 
the scattering rates, rj(x) = 0. All other components 
then remain those of the unmodulated system. If, on the 
other hand, there is only a modulation of the scatter- 
ing rates while V(x) = and B m (x) = 0, only p yy will 
be affected. To obtain (within the classical Boltzmann 
equation approach) a modulation effect on more than one 
component of the resistivity trensor, one needs a mobility 
modulation in addition to an electric and/or a magnetic 
modulation. We will make explicit use of the equation 
of continuity to derive these results, which generalize an 
early observation by Beenakker [|l0| to more general types 
of modulations and scattering mechanisms. Beenakker 
solved the problem for a purely electric modulation in 
the relaxation time approximation. To appreciate the 
important role of the continuity equation in this context, 
it is instructive to consider first the simple local limit. 

We want to emphasize that we assume translational 
invariance in the y-direction throughout this work, so 
that, together with the modulations, current densities 
and electric field components may depend only on x. 
Boundary effects are neglected. The experimental situa- 
tion of a Hall bar with current flow in x-direction will be 
simulated by assuming that the average current density in 
y-direction vanishes. The experimental boundary condi- 
tions of exactly vanishing j y (x) on the sample boundaries 
in y-direction would destroy the translational invariance 
in y-direction and, thus, the simplification resulting from 
the merely unidirectional modulation. They would re- 
quire a fully numerical treatment already in the local 
limit JTgj ] , which may become necessary under certain 
circumstances, but will not be discussed in the present 
work. 



A. Local limit 

In the local limit we describe the linear relation be- 
tween current density j(x) and driving electric field F(x) 
by 

p loc (x)j(x)=F(x), (3.1) 

with a local resistivity tensor p loc (x). It is obtained 
from the Drude resistivity tensor p D of the homoge- 
neous system, with components p xx — p yy — po and 
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p xy = -Pyx = ^cftrpo, where po = m/(e 2 T tr n e i), by 
replacing the spatially constant electron density n e i, cy- 
clotron frequency uj c , and transport time ft r by their spa- 
tially modulated counterparts n e i (x) = n e ;[l — V(x)/Ep], 
w(x) — lo c + u} m (x), and T tr (x) = f tr /[l + f tr r tr (x)], re- 
spectively. In thermal equilibrium, Eq. ( |3.l|) applies with 
j(x) = and the vanishing electrochemical field F(x) 
given by the sum of the electrostatic field W(x)/e and 
the "chemical" contribution Vn e ;(x)/(e-D ). The sta- 
tionary response to a homogeneous external electric field 
E° is accompanied by a change Sn e i(x) of the electron 
density, so that F(x) = E° + VSn e i(x)/(eD ) has to be 
Since Sn e i (x) must have the peri- 



inserted into Eq. (3.1 



odicity of the electric modulation, we find (F(x)) = E°, 
where (...) denotes the average over a period of the mod- 
ulation. Since current density and field depend only on x, 
in the stationary state the continuity equation V ■ j = 
and Maxwell's equation V x F = imply that j x and 
F y must be independent of x. With this result, we can 
define an effective resistivity tensor p by 



p-(j(x))=E° 



(3.2) 



and calculate it by considering two experimentally rele- 
vant situations. First, we assume that the average cur- 
rent density in y-direction is zero. Then, solving the y 
component of Eq. (3.1) for j y (x), we obtain 



Fy 

Pyx — . 

Jx 



I loci lo 
\ryx I ryi 



(3.3) 



yield any magnetoresistance. If only the mobility is mod- 



ulated, one has (1/p 1 . 



loc\ 



Co {Ttr(x))/rt T and only p yy 



deviates from the corresponding Drude value for the un- 
modulated system. If all types of modulation are present, 
all components of the effective resistivity tensor deviate 
from the Drude values. For weak modulations, one ob- 
tains up to second order in the modulation amplitudes 



Pxx Pxx 

Po 



Pxy Pxy 
PO 

Pyy ~ Pyy 



■ v Y ( v 

+ ( ftrnr jk)' 



V 




(3.7) 



"/, ''/, | ^'rT tr — + TtrUJ m ) ) , (3.8) 



Po 



x hip 



(3.9) 



We want to emphasize that, in order to obtain these 
results, it was important to satisfy the equation of conti- 
nuity. If instead we would have taken naively the spatial 
average of the local conductivity tensor, the resulting ef- 
fective resistivity tensor would not have the correct struc- 
ture. A pure electric modulation would have no effect at 
all. A pure magnetic modulation would affect p xx and 
Pyy in the same manner. Although the leading order cor- 
rection to the xx-component would agree with Eq. (ft.7|) , 
the corresponding correction to t he ^ -component would 
disagree with the correct result (3.9). 



Inserting j y (x) into the x-component of Eq. 
obtains 



one 



B. Nonlocal calculation 



(Fx) 

jx 




\o - n loc ] 

IPyX Pyxl i oc 

ryy i 



(3.4) 



In the second, even simpler calculation, we set j x = 
and obtain 



Pyy 



3l 

(Jy) 



1 



and 



Pxy 



(Fx) _ (P l x°y c /P' 

(Jy) 



(3.5) 



(3.6) 



loc 



Since = -p xy 



one finds from Eqs. (3.3) and (3.6) 
Pyx — ~Pxy If there is no mobility modulation, r(x) = 0, 
one has (l/p l yy) — era with cto = 1/po the Drude conduc- 
tivity of the homogeneous 2DES. Since (u>(x))f tr = u> c f tr , 
Eqs. (3.3), fl3.6|), and (3.5) reduce to the correspond- 
ing components of the Drude resistivity for the homog- 
neous system. Only p xx = (p l ™{x)) +?l[(u(x) 2 p l ™(x)) - 
uj 2 /(To] is different from the Drude value l/co, and shows 
a positive magnetoresistance (oc lo 2 ). Note that this re- 
sult differs from the simple average (p l xx ) , which does not 



The local limit applies to systems in which the elec- 
tronic mean free path is much shorter than the modu- 
lation period. In order to exhibit the Weiss oscillations, 
the system must be in the opposite limit, with a mean 
free path considerably larger than the modulation pe riod . 
For such systems, we consider the current density (2.S) 
resulting from the solution of Boltzmann's equa tion (2.7) 
and define the effective resistivity tensor by Eq. ( ft.2| ) . We 
want to emphasize that the spatial average affects only 
j y (x), whereas j x must be independent of x by virtue of 
the continuity equation. 

By formal, but exact manipulations of the Boltzmann 
equation we show in Appendix |X| that the tensor inverse 
of p can be written as 



Pqo 



D f& D 



(3.10) 



where a D is the Drude conductivity tensor of the cor- 
responding homogeneous system, and F has the compo- 
nents 



r, 



with 



-(T tr r tr V/E F ) - G xx , r a 



-(T tr r tr V/Ef 



Gy 



-G x 

'Gy 



(3-11) 
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r - ^IlL 

^ u - 2 
F 



dtp 
2^ 



Here 



(I, 



vp dV/dx 
~2 E v 



+ UJ m v y , g y = —r tr Vy , 



(3.12) 



(3.13) 



and x v is the solution of the modified Boltzmann equa- 
tion 



[V -C] X u = 9u ■ 

Then, by tenso r inversion, a 1 = 
from Eq. ( |3~Tc| ) 



1 + Po a u T 



(3.14) 

p D + A/3, we obtain 



(3.15) 



Explicit matrix inversion and multiplication finally leads 
to the relatively simple result 



Ap 



dctf p D 



l + Cj2 T 2, p Q 



1 



(3.16) 



with detr = T xx T yy - T xy T yx and 



^ _ | _|_ r xx -+- Tyy + u! c Tt r (T xy Tiyx) + detr . ^ 



1+UJ 2 T? T 



is the o nly n onzer o ma trix element of T, 



if r^ y or r,__ 

Eq. ( |3.16 ) reduces to ( |A16| ) or (A17), respectively, i.e., 
either only Ap xx or only Ap yy is nonzero, and we obtain 
the same tensor structure as in the local limit. As shown 
in Appendix [A| these resul ts can be derived without the 
explicit formal result (3.16). 

We want to emphasize th at on ly the global transport 
time fi r , defined below Eq. ( 2.12 ), enters these consider- 
ations about the structure of the effective resistivity ten- 
sor, not any further details of the scattering mecha nisms . 
The latter enter, of course, the s olutions of Eq. ( 3.14 ) 
and determine, according to (3.12), the numerical values 
of the tensor components of T. 



order, we should negle ct m odulation effects on C and on 
V, i.e. replace in Eq. ( |2.9| ) v(x) with vp and omit ui m (%) 
and u) e (x, tp). In this approximation the different Fourier 
coefficients of the expansions V(x) = J2 q9 Lo ^Q ex P{^Q x )i 
and similar for u> m (x) and rt r [x) with coefficients uj q and 
r q , respectively, do not mix in Eq. ( |3.14| ) and lead to un- 
coupled integro-differential equations of the variable (p. 
We now follow Beenakker pof and consider only isotropic 

scatterin g, i.e . we take Jq = 1 and 7™ = for m ^ 
in Eqs. ( |2.11 ), and define the relaxation time r = f by 
the spatial average 1/r = (l/r(x)). Then we can solve 
these equations analytically and obtain 



G *»= ^E {\MV q /E F \ 2 Ql 



(4.1) 



9#0 



- [A\qru q V- q lE F ]Q e ™ + \2ruj q \ 2 Q q nm } , 
G xy = J2 Tr- q [{\qV q /E F )Q q m - 2r^Q™ m ] , (4.2) 



q^O 



Gy X — — G X y , Gyy — 2 llT^C?™™ , (4.3) 

q^O 



with Q^ e = S q m /N q , Q e q m = S^ 1] /N q , and C 
S q hl) + (S q °' 1] ) 2 /N q , where N q = 1 - 4°' 0) and 
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q 



oo 

E 

m— — oo 



J&\Rq)J&\Rq) 
1 + (tocj c t) 2 



(4.4) 



where A = Wi?T is the mean free path, i? = vp/u> c is 

the cyclotron radius, and Jm\z) — J m (z) and Jm^z) = 
J' m (z) are the Bessel functions and their derivatives. For 
uj q = 0, r q = 0, and Vg = 5\ q \ t 2 n / a V r ms/V2 this reduces 

to Beenakker's result [[[o]. The term S q ' ^ in the de- 
nominator N q arises from the back-scattering term in the 
collision operator, which describes relaxation towards the 
local equilibrium. It is thus a consequence of the equation 
of continuity. 



1. Clean limit, uj c t 2> 1 



IV. ANALYTIC RESULTS 

A. Weiss oscillations 

If the modulation is weak, in the sense that the forces 
due to the electric and the magnetic modulations are 
weak as compared to the Lorentz force due to the av- 
erage magnetic field Bo, a perturbation expansion with 
respect to the amplitudes of the modulation is possible. 
To obtain the conductivity correct to second order, we 



need, according to Eq. (3.12), Xv to firs t ord er, and we 
can replace the square bracket in Eq. (3.15) by unity. 
Since the right hand side of Eq. (3.14) is already of first 



In the limit lo c t 1, only the term with m = 
survives in the sum of Eq. (4.4), and in Eq. (A16) the 



term T xx may be neglected, so that Ap xx /po « G xx . 
Then, r earr anging terms with q and —q, the curly bracket 
in Eq. (4.1) can be replaced by \XqS q /Ep\ 2 with S q = 
V q J (Rq) + {k F /q)Tuj q J x {Rq), 



Ap xx /p 



2 ^ 1 



\XqS q /Ep\ 



q^O 



(4.5) 



This result is, apart from the denominator 1 — J 2 (Rq) , in 
agreement with a recent simplified calculation fllq] based 
on an evaluation of Chambers' formula for the conduc- 
tivity in terms of the drift velocity of cyclotron orbits. 
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That calculation corresponds to a naive relaxation time 
approximation in Boltzmann's equation with relaxation 
towards the total instead of the local equilibrium distri- 
bution function, and thus violates the equation of conti- 
nuity for the modulated system. Indeed, the nontrivial 
denominator in the present calcu lation results from the 
back-scattering term in Eq. ( [2.10 ). 

In Fig. [I] we demonstrate the range of validity of sev- 
eral approximations to the results (4.1) - (4.4) for the 
special case of a simply harmonic electrical modulation 
V(x) — Vocosqx. Figure [l](a) shows, for several values 



of the mean free path A, the quantity Sq ' ^ of Eq. (4.4), 
which determines the denominator of the magnetoresis- 
tance corrections. Figures 0(b) and (c) show resu lts for 
the magnetoresistance, calculated from Eq. ( |4.l| ), with 
different normalizations which are suitable to discuss the 
local limit, (b), and the limit of infinite A, (c), respec- 
tively. The curves are plotted versus the inverse cy- 
clotron radius, which is proportional to the average mag- 



netic field Bn. In the limit A — > oo, S, 



(o,o) _ T 2 



(2 / ir Rq) cos 2 (Rq — 7r/4) decreases with decreasing B , 
with an infinite number of zeroes and with values at the 
maxima, which approach zero linearly with Bo. The same 
then hold s for the magnetoresistivity curve obtained from 
Eq. (fh^) and for the thin line in Fig. [j] (c), which cor- 
responds to the simp lifie d result of Refs. JT^jT^] and is 
calculated from Eq. ( |4.5| ) by replacing the denominator 
by 1. For 1/ 1 Rq < 0.1 the denominator of Eq. (4.5) is close 



to unity (S^ 0,0 ' < 0.1 ), and the thin line approximates 
the result of Eq. (4.5) quite well. At higher magnetic 
fields, l/Rq > 0.2, the result of Eq. ( fOj| ) agrees well 
with the magnetoresistivity curve obtained for Xq = 30 
in Fig. |] (c). For a sam ple with a given value of the 



mean free path, Eq. (4.5) does not hold uniformly for 



all values of the average magnetic field. With decreasing 
Bo values, Xq/Rq — ui c t becomes smaller an d more and 
more terms contribute to the sum in Eq. (4.4). As a con- 
sequence, the minima of the Weiss oscillations of S q ' ^ 
and Ap xx are risen to positive values, the amplitudes of 
the oscillations are increasingly suppressed with decreas- 
ing Bq values, and, in the limit Bq — > 0, a finite positive 
value is approached with zero slope. 

Whereas Sq ' ^ in the limit of large Bq (Rq — ► 0) 
approaches 1 for all values of A [since Jq(0) = 1 and 
Jm(0) = for m 7^ 0], it shows an overall increase 
with decreasing A for all values of Rq, and approaches 
1 uniformly in Bo for Xq — > [since Y]^^ Jm ( z ) = !]■ 

in Eqs. (p 



Thus, the denominator N„ 



-(4.3) becomes 
most important in the local limit Xq — > and in the 
"quasi-local" limit Rq < 1, where the cyclotron radius 
R becomes smaller than the period a = 2n/q of the 
modulation. Figure [j] (b) shows the magnetoresistivity 
data for Xq = 9 and Xq = 4 in a representation which 
is more convenient for the discussion of the local limit 
Xq — » 0. For Xq = 1, the corresponding result ap- 
proximates already closely the local result of Eq. (3.4), 
Ap xx /p Q = {\/2){Vo/E F ) 2 [l + (X/R) 2 }, which is indi- 



cated by the dash-dotted line. 

In the "quasi- local" high magnetic field limit, Rq 
1, the denominator becomes 1 — Jo(Rq) « (Rq) 2 /2 <C 
1, so that, for V(x) = Vbcosgx, we obtain a quadratic 
magnetoresistance, 

Ap xx /p « {u c r) 2 (Vo/E F ) 2 /2 , {Rq « 1) , (4.6) 

whereas the calculation of Ref. Jlfj] yields a saturation 

for large uj c , Ap x l x / Po ~ {M) 2 (Vo/Ef) 2 /^. Experiments 
seem to show a quadratic magnetoresistance at high mag- 
netic fields, although masked by the onset of Shubnikov- 
de Haas oscillations. For a sinusoidal magnetic modu- 
lation, uj m {x) — cj°„ cosqx, we obtain from Eq. (15) at 
large Bo values a saturation, 



Ap xx /p Q ^(uj m T) 2 /2, (Rq<l) 



(4.7) 



whereas Ref. @ yields A P p/p w (A<z) 2 «/w c ) 2 /4, 
which becomes small oc Bq 2 . We want to emphasize 
that, in this "quasi-local" limit, the Bq dependence of the 
magnetoresistance obtained fro m E q. ( p~5| ) agrees with 
that of the local approximation (3.7). 

Figure || shows the prediction of Eqs. ( [4.l|) - (4.4) for 



the case when two types of modulations are simultane- 
ously present. In Fig. || (a) an electric and a magnetic 
cosine modulation of comparable effective strengths are 
considered. The dash-dotted line refers to a pure elec- 
tric and the solid line to a pure magnetic modulation, 
whereas the dashed lines refer to a combination of both. 
For fixed amplitudes of the electric and the magnetic co- 
sine modulations, the resulting magnetoresistivity curve 
depends strongly on the phase difference between both. 
The long-dashed line is Fig. [2] (a) is obtained for in-phase 
modulations, and the short-dashed line for a phase shift 
of a half period. A systematic investigation of the in- 
terference effects occurring due to the superpositon of an 
electric and a magnetic modulation has been given within 
the simplified treatment of Ref. Jl5| . Since for Xq > 20 
this treatment is qualitatively correct in the regime of 
Weiss oscillations, we will not repeat this discussion here. 

Figure ^| (b) shows the components of the magnetore- 
sistivity tensor for a combination of electric and mobility 
modulations. We assume here in-phase cosine modula- 
tions for the electrostatic potential energy and the scat- 
tering rate of the electrons, and c omm en t on opposite 
phases below. According to Eqs. (4.1) - (4.4), Ap xx is 



determined by the quantity Q q e , Ap yy by Q r q m , and 
Ap xy — —Ap yx by Q q m , all of which contain different 
combinations of Bessel functions and their derivatives. 
For Xq > 20, in the regime of Weiss oscillations Q q e is pro- 
portional to the square of Jo(Rq) and Q q " n , which also 
determines the results of a pure magnetic modulation, 
is proportional to the square of the derivative J' (Rq). 
This leads to antiphase oscillations of Q q e and Q™ n as 
depicted in Fig. || (a) by the dash-dotted and th e soli d 
lines, respectively. However, according to Eqs. ( 3.11 ), 



G 



( |3.15 ), and (4.1) - (4.3), Ap xx and Ap yy exhibit in- 
phase oscillations since, apart from i?o-independent off- 
sets +rr° n (Vo/E F )/2 and -Tr° n (V /E F )/2, respectively, 
Ap xx is proportional to Q 3e , whereas Ap yy is propor- 
tional to — Q™ m . Thus, a mobility modulation can not 
explain the anti-phase oscillations of Ap yy observed in 
the early experiments [[l], as mentioned in the Introduc- 
tion. In the situation considered in Fig. |[ Q 3m is ap- 
proximately given by the derivative of Q 3e with respect 
to Rq. Thus in Fig. ^ (b), where the resistivity com- 
ponents are plotted versus 1/Rq, Ap yx appears as the 
derivative of Ap xx . If we change the relative sign of elec- 
tric and mobility modulation, the sign of Ap yx will also 
change, whereas the phase relation between Ap xx and 
Ap yy remains unchanged. 

We do not show a figure analog to Fig. || (b) for a 
superposition of a magnetic and a mobility modulation, 
since in this case all components of the magnetoresistivity 
tensor are determined by Q™ m . As a consequence, Ap yx 
would show the same appearance as Ap XXl i.e. as the 
solid line in Fig. || (a), and Ap yy would appear as —Ap xxi 
i.e. with the opposite phase. In contrast to the case of 
electric and mobility modulation, in the case of magnetic 
and mobility modulation there is no constant offset in the 
diagonal components of the resistivity tensor, and Ap yx 
does not assume negative values. 



2. Damping of Weiss oscillations 

The value of Sq ' ^ for weak magnetic fields, Rq > 
2, depends strongly on Xq, and Weiss oscillations are 
washed out for 1/Rq < 0.5/Xq, i.e., for u c t < 0.5. If 
the mean free path is too small, Xq < 2 (here we have 
in mind that the period of the modulations, and thus q 
is fixed), no oscillations survive, since too many Bessel 
functions contribute to the sum in Eq. (4.4) (ui c t < 1 in 
the 1 / Rq regime where oscillations can be expected) . In 
this regime the calculation of Ref. [|l5| is not applicable, 
since it assumes X ^> R. 

We can al so r e-obtain the local limit Xq — > from 
Eqs. (4.1) - (O), if we expand the Bessel functions in 
Eq. (4.4) for small arguments, Rq — > 0, and keep only the 
leading orders. This means that we consider the limit of 
large modulation period, but make no assumption about 
the magnitude of uj c t. In this limit, only the terms with 
m = and m = ± 1 co ntr ibut e to Eq. (4.4), and we ex- 
actly recover Eqs. ([Tt]) - (3J)), with t instead of f tr . To 
obtain this result correctly, it is essential to calculate the 
denominator (which vanishes in this limit) consistently 
up to the order (Rq) 2 , which includes contributions from 
the Bessel functions Jo an d J\ ■ Th is d emo nstrates nicely 
that the denominator in Eqs. (LI) - (O) indeed is im- 
porta nt in or der to obtain the correct local-limit results, 
Eqs. ( |3.7| ) - (3J3), which were derived by exploiting ex- 
plicitly the equation of continuity. 



B. Low- field magnetoresistance 

If the average magnetic fi eld B p becomes too small, 
the analytic solution of Eq. (3.14) for small modulation 
amplitudes is no longer applicable. It was obtained un- 
der the assumption, that for the calculation of \v the 
modulation contributions u) m (x) and u> e (x, tp) to the drift 
operator ( |2.9P could be neglected as compared to the av- 
erage cyclotron frequency ui c . Physically this means that 
the modulations modify the cyclotron motion in the av- 
erage magnetic field only weakly and lead essentially to 
a drift of the cyclotron orbits. If however, for a given 
modulation, Bq becomes smaller than a critical B cr a (or 
R > Rcrit), this is no longer true and other types of or- 
bits ( "channeled orbits" ) occur which are confined in the 
x-direction to a single period of the modulation super- 
lattice. The effect of these orbits, which lead to a posi- 
tive magnetoresistance at small Bq [|J , is not included in 
Eqs. (Q - (fOI). 



1. Classification of orbits 

To characterize the channeled orbits and to elucidate 
their importance for the magnetoresistance, we proceed 
as follows. We note that, for a one-dimensional mod- 
ulation in ^-direction, there exists a se cond conserved 
quantity in addition to the energy, Eq. (^^.Integrat- 



ing the y-component of Newton's equation (2.1), we can 
write 



e 

mc 



xBn 



dx'B m (x') 



Py 
rn 



(4.8) 



where p y is a constant of motion, which has the meaning 
of a canonical mo mentum. Solving for v y and inserting 
this into Eq. (^^), we obtain for the projection of the 
orbit on the a;-direction 

(m/2)v 2 x + V(x) = E F , (4.9) 

with the effective one-dimensional potential 

,B m (x>y 2 



V(x) = V(x) + -u? c 



X — Xq + I dx 

B o 



(4.10) 



Here we assume Bq ^ 0, and the constant of motion is 
written as Xq = —cp y /eB , which in the absence of mod- 
ulations is the x coordinate of the center of the cyclotron 
orbit. For each value of xq one obtains a different location 
of the effective potential V(x) on the x-axis, and thus a 
different orbit located between turning points, at which 
v x = 0. If a; is a turning point, V^a;) = E F , Eq. (4.1C) 
yields two possible values for the constant xq: 



Xq(x) 



dx'^l±R 

Bq 



E F 



(4.11) 
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Since V(x) < Ep holds for any orbit passing through a 
given position x, only orbits with xq values in the inter- 
val Xq (x) < x < Xq (x) can pass through x. (We specify 
orbits only up to an arbitrary shift in y direction.) On 
the other hand, a given value of xq determines an or- 
bit between turning points Xt satisfying x^{xt) — xq or 
x~Q(xt) — xq. Thus the curves xq = x^(x) yield the con- 
stant of motion of the orbits which have i as a turning 
point. We use this property to distinguish two types of 
orbits. 

Orbits of the first type we call "drifting orbits" . These 
have their left turning point xi satifying x q (xi ) = xq and 
their right turning point x r satifying Xq(x t ) = xq. For 
sufficiently smooth modulation (cf. Appendix B) and 
large Bq this is the only type of orbits, as we demon- 
strate in the left panels of Fig. |^ for an electric cosine 
modulation. In the limit of vanishing modulation these 
"drifting orbits" reduce to cyclotron orbits, i.e., to cir- 
cles, and x r — xi — > 2R. The effect of the modulation in 
x-direction is to deform the circular cyclotron orbits and 
to superimpose on the cyclotron motion a drift velocity 
in y-direction. The drift vecolity depends on the position 
of the orbit in the modulation field and may vanish for 
highly symmetric positions. In the left panel of Fig. |3| 
one orbit is nearly symmetric with respect to the poten- 
tial maximum at x = 0, another one nearly symmetric 
about the minimum at x = 1.5a. The small distance 
between neighboring loops of the orbits (three loops are 
shown) indicates a low drift velocity. 

Orbits of the second type, which we call "channeled or- 
bits", occur if the curves Xq (x) have local extrema, i.e., 
for sufficiently strong modulation or weak magnetic field 
Bq. For "channeled orbits" we require that their turn- 
ing points both satisfy either Xq(xi) = Xq(x t ) = xq or 
Xq(xi) = Xq(x t ) = xq, i.e. in plots like those in the up- 
per panels of Fig. || both turning points are located either 
on Xq{x) or on Xq(x). These channeled orbits (see mid- 
dle and right panels of Fig. ||) are confined to a single pe- 
riod a of the modulation and have no counterparts in the 
unmodulated system. For an electric modulation they oc- 
cur near potential minima. For a magnetic modulation 
they occur where the total magnetic field Bq + B m (x) 
changes sign |2l|| . 

Each value Xq of the constant of motion determines 
uniquely a single drifting orbit (of course, modulo a shift 
in y-direction). Channeled orbits may occur with the 
same xq value, however only in intervals of the x axis 
which are not entered by that drifting orbit. At a given 
position x, each allowed xq value characterizes uniquely a 
single orbit through x. In the limit -Bo = 0, drifting orbits 
degenerate into wavy trajectories with periodic velocity 
v(x + a) = v(x) and nonzero mean velocity (v x ) in in- 
direction. The mean value in y-direction is (v y ) = p y /m 
(note that Xq/R = —p y /mvF, and for pure electric mod- 
ulation v y — Py/m is constant). In this limit two trajec- 



If xq coincides with the value of a local maximum of 
Xq(x) or a local minimum of Xq(x), the orbits near the 
corresponding turning point Xt become critical in the 
sense that they do not reach this turning point, but ap- 
proach in the x — y plane the straight line x = x t asymp- 
totically. 



2. Effect of channeled orbits 

As has been emphasized by Beton et al. [Q, the chan- 
neled orbits are responsible for the strong positive magne- 
toresistance of modulated systems in weak perpendicular 
magnetic fields. In Appendix ^ we calculate the fraction 
f co {Bo, Vq) of the phase space covered by channeled or- 
bits, as a function of the (average) magnetic field Bq and 
the modulation strengths Vq, for two electrostatic mod- 
ulation models. For the smooth cosine modulation with 
a fixed Vq, channeled orbits exist only at sufficiently low 
Bq < B cr it(Vo), whereas for a discontinuous step modu- 
lation skipping orbits survive at arbitrarily high magnetic 
fields. For both models the fraction of channeled orbits 
at zero magnetic field, / co (0,Vb), increases at low mod- 
ulation strengths proportional to \Vq/ Ep^^ 2 ■ We also 
sketch, for Bq = and in the nonlocal limit of large 
mean free path, a calculation in the spirit of Ref. |lq | 
which shows that the resistivity correction increases as 
\Vq/Ef\ 3 / 2 with the modulation amplitude. This indi- 
cates that, for Bq = and weak modulation, no expan- 
sion of the transport coefficients in a power series of Vq 
is possible. The same is true at finite magnetic field if 
channeled orbits exist. Indeed, for the step model which 
allows for channeled orbits at arbitrary strong Bq, the 
sum over the Fourier coefficients in the second-order for- 



tories with 



> and 



< occur with the same 



as indicated by the dashed lines in the lower right 



panel of Fig. 



mula (4.1) for the magnetoconductivity diverges for all 
values of Bq, indicating that such an expansion does not 
exist Q. 

To include the effect of channeled orbits on the magne- 
toresistance properly, we have solved Boltzmann's equa- 
tion numerically. The results show in the limit of 
long mean free path and zero magnetic field the same 
\Vq/Ef\ 3 / 2 behavior as obtained analytically in the sim- 
plified approach, and will be presented in the next sec- 
tion. 



V. NUMERICAL RESULTS 

A. Some technical details 

In order to cover the whole magnetic field range, we 
have solved Eq. (3.14) numerically, using Fourier expan- 
sions fo r the dependences on both x and (p . Ac cording 
to Eq. ( 2.11 ), the collision operator, Eq. ( 2.10 ), is di- 
agon al in the Fourier representation. The drift term, 
Eq. (2.9), becomes simple for strictly harmonic modu- 
lations, provided that the electric modulation is weak 
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enough and v(x) can be approximated by vf- Then 
the drift term c ouple s only neighboring Fourier coeffi- 
cients, and Eq. (3.14) leads to an infinite set of linear 
equations for the Fourier coefficients with a sparsely oc- 
cupied matrix. Truncating the infinite set by consider- 
ing only the lowest N x and N v Fourier coefficients for 
the x- and the ^-dependence, respectively, we could ex- 
ploit the sparse nature of the system matrix and handle 
approximations of the infinite set with large dimensions 
N x ■ N v . The number of Fourier coefficients needed to 
obtain converged results increases with increasing mod- 
ulation strenghts and with increasing mean free path. 
The convergence in the regime of Weiss oscillations is 
much better than in the low-£?o regime, where the chan- 
neled orbits are important. For the pure electric cosine 
modulation with Vq/Ef = 0.02 we needed N x = 13 and 
N v = 550 to achieve convergence in the low-magnetic- 
field regime 5 ■ 10~ 5 < 1/qR < 0.02 , whereas N x = 4 
and N v = 450 was sufficient for 1/qR > 0.5. The reason 
for the poor convergence at low B is th at the distribu- 
tion function calculated from Eq. ( [3.14 ) developes very 
sharp and rapidly changing structures near the angles 
ip = 7r/2 and 37r/2, originating from the turning points 
of the channeled orbits. Owing to these numerical prob- 
lems, our Fourier expansion approach is reliable only for 
weak modulations (Vq/Ef < 0.2). 



B. Relaxation time approximation 

A typical low-magnetic-field result for the pure elec- 
tric modulation V(x) — Vbcosqx with Vq/Ef = 0.02 is 
shown in Fig. |i| The solid line is the result of our numer- 
ical calculation. The dotted line shows the result of the 



analytical approximation (4.1). The dash-dotted line in- 
dicates the contribution of the channeled orbits and is ob- 
tained as follows Q. For weak modulation and lo c t S> 1, 
the correction to the Drude resistivity tensor can be writ- 
ten as Ap xx /p w uj^T 2 A(Tyy/(To- The contribution of 
channeled orbits may be estimated by their m ean s quare 
drift velocity Aa yy /a = (vf)/(v 2 F /2), cf. Eq. flBl^) . Ap- 



proximating v d s» v F , this becomes » 2/™ s (_B , Vq), i.e., 
twice the fraction of channeled orbits at the Fermi energy 
(see appendix |bJ). The dash-dotted line in Fig. | repre- 
sents Ap xx /p w 2uj 2 T 2 fc° s {B ,Vo) with an additional 
vertical shift to match our numerical result at Bq = 0. 
The vertical straight line indicates the critical magnetic 
field Bf r jj° n = cqVo/evp at which the channeled orbits 
die out . Our numerical result exhibits a maximum at 
the magnetic field value B p sa 0.8B^ r f t on . This is in very 
good agreement with recent experiments p0| , although 
the numerical factor 0.8 should depend on the specific 
form of the potential modulation. The corresponding re- 
sults for a pure magnetic modulation look qualitatively 
very similar and are not shown. 

The inset of Fig. ^ shows our numerical results for 
strong purely magnetic modulations of the form B m (x) = 



B^ cos qx. The modulation strength is given in units of 
1/qRm = aeB { ^ n /{2TTmcv F ), Ap xx /p in units of w^r 2 = 
X 2 /R 2 l . With increasing modulation strength (different 
curves in the inset of Fig. [|) , the extent of the positive 
magnetoresistance regime at low Bq values increases and 
Ap xx reaches a maximum at B p m 0.7 B^ (note that 
B C rit = B^, see appendix |b|). The magnetoresistiv- 
ity at Bp increases with increasing modulation strength, 
and the low-field Weiss oscillations are successively sup- 
pressed (the fluctuations in the curve for l/qR m = 0.34 
are due to numerical errors, which become larger with 
increasing modulation strength). Again the dependence 
of the magnetoresistance on the modulation strength is 
very similar for the electric cosine modulation and will 
not be shown explicitly. For large mean free path and 
Vq/Ef > 0.2, fluctuations of the calculated curves indi- 
cate convergence problems. We found, however, that the 
relation B p ps 0.8B^l° n , i.e. a linear increase of B p with 
the modulation amplitude Vq, is reasonably well satisfied 
up to Vq/Ef ~ 0.5, in agreement with Ref. M and with 
a recent modification of Stfeda's semiclassical approach 
fL3|. The peak values of the magnetoresistivity can be 
well fitted by 

Ap xx (B p )/p Q = 7.01 (Vq/Ef) 2 + 54.7 (Vq/Ef? . (5.1) 

Whereas in the regime of Weiss oscillations (Bq > 3B cr n) 
Ap xx is prop orti onal to (Vq/Ef) 2 and the analytic ap- 
proximation (4.1) is valid, the peak va lue at B p is domi- 
nated by the third order term in Eq. (5.1) for Vq/Ef > 
0.15. For weak modulation (Vq/Ef < 0.15), Xq^> 1 and 
Bq = 0, our numerical results are consistent with the an- 
alytically calculated algebraic dependence Ap xx (0) / p cx 
(Vq/Ef) 3 / 2 (see Appendix |b|). This indicates that, at low 
magnetic fields where channeled orbits exist, an expan- 
sion of the magnetoresistance in powers of the modula- 
tion strength is not possib le. T o check this, we expanded 
the Boltzmann equation ( [3.14 ) in powers of Vq/Ef and 
solved it up to the order (Vq/Ef) 8 - In the regime around 
B p the higher order contributions lead to high-frequency 
oscillations, but we found no indication of convergence 
or positive magnetoresistance. 



C. Anisotropic scattering 

Experiments in the regime of Weiss oscillations give 
evidence that the mobility of the 2DES is dominated by 
scattering of the electrons by remote ionized impurities, 
which should lead to predominant forward scattering. To 
investigate the effect of anisotropic scattering, we simu- 
late the angular dependence of the different ial sc attering 
cross section in the collision operator, Eq. ( |2.10 ), by the 
simple model 



P p (<p) = c+(l-c) 



n 2 



(2p)\ 



2cos-\ 
2) 



2p 



(5.2) 



9 



where c controls a constant background (0 < c < 1), 
and p the sharpness of the forward scattering peak. This 
model is convenient, since it has simple Fourier coeffi- 
cients, 70 = 1, 7,„ = (1 - c)(p\) 2 /\(p + m)!( p - m )!] for 



\m\ < p, and 7™ =0 for |m| > p [cf. Eq. ( |2.1l| )]. Ac- 
cording to Eq. ( [2.12 ), the ratio of transport time and 
relaxation time r increases then with increasing p as 

T tr /T=(p+1)/{C P +1). 

Keeping r fixed, we have calculated magnetoresistiv- 
ity curves for several values of the parameters c and p. 
We found as a general trend that, with inreasing r tr ./r, 
for both electric and magnetic modulations the values of 
Ap xx /po [with po = m/ (e 2 n e iTtr)] increase slightly faster 
than proportional to t 4i ./t. On an absolute scale this 
means that the values of Ap xx / pq change little near well 
developed minima but increase stronly at the peaks, re- 
sulting in a strong increase of the amplitudes of the Weiss 
oscillations. In Fig. |5| we have, for increasing T tr , re- 
duced the relaxation time r so that the maximum value 
of Ap xx /po near 1/qR = 0.26 remains unchanged. In 
this plot higher-index Weiss oscillations are increasingly 
damped with increasing importance of small-angle scat- 
tering. Comparison with typical experiments shows that 
small-angle scattering is indeed important. Apart from 
the quantum mechanical Shubnikov-de Haas oscillations 
at large Bq, we could nicely fit experimental curves in the 
whole region 1/qR < 0.8 with typical values of r tr /r > 2. 



D. Mobility modulation 

Finally, we have considered, within the relaxation time 
approximation, a mixed electric modulation V(x) = 
Vocosgx and mobility modulation 1/t(x) = f~ 1 + 
r m cos qx. We choose an in-phase modulations, Vq • r m > 
0. We found that, for Vq/E f <C 1 and fVFq ^S> 1, the 
effect of the mobility modulation on p xx is rather small. 
This is understandable in the regime of Weiss oscillations, 
since, according to Fig. ||, Ap xx /p becomes independent 
of Xq for Xq > 25. It is also expected for Bq — in 
the local limit, where Ap xx /po becomes proportional to 
(Vo/Ep ) • fr m (For the corresponding step modulations 
one has p xx /p = [1 + (V Q /E F )fr m ]/[l - {V /E F ) 2 }.) 

In Fig. ^| we show the results for pure electric and pure 
mobility modulations as solid lines. Whereas the elec- 
tric modulation enhances p xx , the mobility modulation 
reduces p yy . The Weiss oscillations of both resistivity 
components are in phase, as w e have seen already from 
the analytical solution Eq. (4.1) - (4.3). For the pure mo- 
bility modulation, there is no positive magnetoresistance 
at small Bq values, since there are no channeled orbits in 
this case. 

Results for the case of mixed modulations are shown as 
broken lines. For the diagonal resistivities, the changes 
compared to the pure cases are essentially a positive off- 
set of p xx and a negative offset of p yy , as expected from 
Eqs. (3.11) and (3.16). The enhancement of p xx at small 



Bq values (R > R cr it) is plausible, since now in the 
regime of channeled trajectories the mobility is increased. 
Therefore, the importance of these trajectories and thus 
the value of the magnetoresistivity is enhanced. 

For the mixed case, we obtain also a correction to the 
Hall resistivity. To accommodate this correction in the 
same figure, we have subtracted the Drude value, which 
increases linearly with Bq. In the regime of Weiss os- 
cillations, Ap xy is roughly proportional to the slope of 
Pxx/ Po as function of Bq oc 1/qR. 



VI. SUMMARY 

Using the linearized stationary Boltzmann equation 
and including anisotropic scattering in the collision oper- 
ator, we have investigated the effect of one-dimensional 
periodic modulations on the resistivity tensor of a two- 
dimensional electron system in a perpendicular magnetic 
field Bq. Even in the relaxation time approximation 
(for isotropic scattering) we have carefully considered the 
back-scattering term of the collision operator, which en- 
sures relaxation of the distribution function towards its 
local equilibrium value in the stationary state, and not 
to the total equilibrium distribution. This is of principle 
importance for the modulated systems, since it ensures 
that the approximation does not violate the equation of 
continuity, and it has qualitative as well as quantitative 
consequences for the resistivity tensor. 

The qualitative consequences concern the form of the 
resistivity tensor. In the absence of mobility modula- 
tions, electric and magnetic modulations in ^-direction 
affect only p xx , whereas the other components of the re- 
sistivity tensor remain those of the unmodulated system. 
This result, which is not obvious in the presence of a mag- 
netic field, has first been obtained by Beenakker |h]] for a 
pure electric modulation in the relaxation time approx- 
imation. In Appendix A we show that this result also 
holds for mixed electric and magnetic modulations and 
for anisotropic scattering. Moreover, we show that, in the 
absence of electric and magnetic modulations, a mo bility 
modulation in ^-direction affects only p yy . In Sect. [II A 
we have derived the same tensor structure in the local 
limit. This more transparent derivation elucidates the 
important role played by the continuity equation. 

Simpler calculations based on an evaluation of the fa- 
miliar velocity-velocity-correlation formula for the con- 
ductivity tensor are not in agreement with these exact 
results on the structure of the resistivity tensor |I^,|l^.|l5| . 
although they may yield reasonable approximate results 
under certain conditions. 

For weak modulations, isotropic scattering, and suffi- 
ciently high magnetic field Bq, we followed an approxi- 
mation scheme proposed by Beenakker |Io|| and obtained 
analy tic results for the magnetoresistivity tensor (Sect. 
IV A ). In the clean limit (mean free path much larger 
than modulation period) these results reduce to those 
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of the simpler calculations, apart from a characteristic 
denominator, which can be traced back to result from 
the back-scattering term in the collision operator. The 
absence of this denominator in the simplified treatment 
Jl4],[l5| indicates that this treatment does not correctly 
describe the relaxation towards the local equilibrium, and 
thus violates the equation of continuity. Quantitatively, 
the deviation of the denominator from unity is not large 
if the cyclotron radius R (of an electron at the Fermi 
energy) is larger than the period a of the modulation. 
Therefore, in the regime of Weiss oscillations the simpli- 
fied treatment gives reasonable results. However, in the 
quasi-local regime of high magnetic fields (R < a), the 
denominator becomes small and the simplified treatment 
becomes wrong. In this regime the correct treatment of 
the equa tion of continuity is quantitatively important, 
see Sect. 



[VA1 



The analytic approximation, which implicitly assumes 
weakly perturbed cyclotron motion, breaks down at low 
magnetic fields, where "channeled orbits" become impor- 
tant j|. At B = and for Xq » 1, these should lead 
to an interesting algebraic dependence of the resistivity 
on the modulation strength, as we demonstrate in Ap- 
pendix B. To include the effect of the channeled orbits, 
we have solved the linearized Boltzmann equation nu- 
merically. The numerical results in the relaxation time 
approximation (isotropic scattering) give a qualitatively 
correct description of the experimental result for p xx in 
the range of small magnetic fields, where a pronounced 
positive magnetoresistance is obtained, and of interme- 
diate magnetic fields, where the Weiss oscillations are 
observed. At Bq = 0, they are also consistent with ana- 
lytical results which are derived in Appendix B. A nearly 
quantitative description of the experimental p xx 0] in 
the weak and intermediate magnetic field regime is ob- 
tained from a numerical calcula tion considering predom- 
inant forward scattering (Sect. V C ). 

In an attempt to find a classical explanation of the 
anti-phase Weiss oscillations observed in p yy [[l] , we have 
also investigated a mobility modulation. We obtained, 
however, only weak in-phase oscillations, both analyti- 
cally and numerically. Thus, this attempt of a classi- 
cal explanation of the p yy oscillations fails, whereas the 
quantum treatment provides a convincing explanation of 
these oscillations as a density-of-states effect 0. 
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APPENDIX A: MOMENTS OF BOLTZMANN'S 
EQUATION 

To investigate the st ruct ure of the tensor p, we multi- 
ply the linearized Eq. (2.7) with e 2 D v(x)u(p) and av- 
erage with respect to both ip and x over a full period. 
Integrating by parts eliminates the der ivatives of the dis- 
tribution function <j)(x, p), see Eq. (2.9), and the collision 
operator is evaluated using the real and imaginary part 
of 



1 



T tr (x) 



* ^e^(x,p). 



(Al) 



where only the effective transport scattering rate l/rt r (x) 
enters [cf. Eq. (2.11)]. Multiplying with f tr and rearrang- 
ing terms we can write the result as 



p D (j(x))=-E° + A, 



(A2) 



where p D is the Drude resistivity tensor of the homoge- 
neous system, and [with V'(x) = dV/dx] 

A* = -{[V\x)/E F ]4> loc {x) + poftrV m (x)j v (x)) , (A3) 



-PQTtr (r tr (x)j y (x)) 



(A4) 



Here we have explicitly used that, according to the equa- 
tion of continuity, j x {x) = (j x (x)) is spatially constant, 
and that the spatial average values of u m (x), V(x), and 
Ttr(x) vanish. 

For the unmodulated system A = 0, and Eq. ( |A2| ) 
reduces to the well known Drude form of Ohm's law with 
homogeneous current density j(x) = = j°. For 

the modulated system with the same E , we may write 
(j{x)) = f + {Sj(x)), so that 



p D (SKx)) 



(A5) 



To calculate these corr ections, it is convenient to write 
the solution of Eq. ( |2.7| ) in the form 



(x, tp) = <j)°{p)v{x)/v F + p T tr x(x, p) . 



(A6) 



where <p (<p) = poTt r vpu(p) ■ j is the corresponding so- 
lution for the unmodulated system. This yields 



with 



j(x)=j°v 2 (x)/v%+6i(x) : 



$j(x) = \ [ v {x)u{p)x{x, p) , 



J F J— it 



2tt 



(A7) 



(A8) 



and x( x i <P) satisfies Eq. ( p/7| ) with a modified inho- 
mogeneity which is linear in the modulation quantities 
w m (x), V(x), and r tr (x), 

[D - C] x = ( - vdv/dx + uj m v y - r tr v x )j x 

- {u> m v x + r tr v y )f , (A9) 
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where {v x ,v y ) = v(x)(cos p, simp). We may write x = 
Xi + X2, where \ 2 solves 

[T> - C] X2 = -«(a0 cos 99 [r tr j° + u m j°] . (A10) 

The solution is independent of tp, and satisfies 

d X 2/dx = ~[rtrf x + Umjy] ■ (All) 

It does not contribute to the current density, Eq. (|A8|) , 
but it contributes to the local equilibrium distrib utio n, 
and thus contributes to the inhomogeneity of Eq. (A3), 



/ dV/dx 
\ E F 



= pon, 



I dV I dx f 7 ' dp 



E F 



2tt 



Xi 



V 



[rtrf x + "mJy. 



(A12) 



for purely electrostatic modulation within the relaxation 
time approach, is thus shown to hold also in the presence 
of a magnetic modulation and anisotropic scattering. 

The other simple case is that without electric and mag- 
netic modulation, V(x) = 0, w m (x) = 0, but with a fi- 
nite mobility modulation r tr (x). In this case A x = and 
A y = poTyyjy. The same procedure yields now a resis- 
tivity tensor, which differs from that of the unmodulated 
system only in the yy component, 

^Pyy/Po = - W(l + T yy /[1 + (iJ c f tr ) 2 }) . (A17) 

In the general case, with mobility modulation and elec - 
t ric or mag netic modulation, we obtain from Eqs. ( (3.2j ), 
(|A2|), a nd (A15 ) the effective conductivity tensor given 
by Eq. fl3.10|). 



From Eqs. ( \Kt\ ) to ( |Al2j ) we obtain 

V ,n 



A x = - /'or,, [ r tr^J 



2 



dp ( vp dV/dx 
E F 



2tt 



+ UmVy Xl 



v .n 2 r dp 

A y = - p T tr { - r tr— J y + — I — r tr VyXl 



E F ~ 



(A13) 



(A14) 



Final ly, w e wri te x.i — Xxj® + Xyi y i where, accordi ng to 
Eqs. (|A9|) and ( |A10| ), Xv (for v = x,y) satisfies Eq. ( §jj ) 
with g v given by Eq. ( |3. 13| ) ■ We obtain 



(A15) 



where f is defined by Eq. ( 3.11 ). 

From these results we can show that the magnetoresis- 
tivity tensor reveals a very simple structure in the two im- 
portant special cases, where there is either no modulation 
of the scattering rate, or where only the scattering rate is 
modulated. Let us first assume that we have an electri- 
cal and a magnetic modulation, but no mobility modu- 
lation, r(x) = 0. Then T xx = — G xx and all other tensor 
components of T vanish, so that = poT xx j®, while 

Ucftr(5j x ) 



A y = 0. Thus, according to Eq. (A5), (Sj y 
and (Sj x ) = T xx j°, with f xx = T xx /[l + (uj c f tr ) 2 }. If we 
now consider the effective conductivity tensor a = p^ 1 
of the modulated system, j° + (Sj) = ctE , we obtain 
for the components of a in terms of those of the Drudc 



v D = [P°] 1; °nv =< T nu( 1 + T xx) for {pv) = (xx), (xy), 



and (yx), whereas a yy = Oy y [l - {u c T tr ) 2 T xx \. Tensor 
inversion shows that p differs from the Drude resistiv- 
ity tensor p D of the unmodulated system only in the xx 
component, 



Ap xx /p = p xx /p - 1 = -T xx /[l + T xx ] 



(A16) 



while all other components are not modified by the mod- 
ulation. This result, first pointed out by Beenakker fl(i| l 



APPENDIX B: CHANNELED ORBITS 
1. Critical field B 

cr 1 1 

For a smooth modulation, the condition for the exis- 
tence of channeled orbits is apparently, that the functions 
Xq(x) and Xq(x) have local extrema, i.e., that the equa- 
tion 



B m {x) R V'{x)/E F 
B T 2 ^1 -V(x)/E t 







(Bl) 



has a solution. For a pure magnetic modulation of the 
form B m (x) = B® n cosqx, we obviously find B c nt = B^. 
For B < B crit the total magnetic field B(x) = B + 
B m (x) changes sign, and channeled orbits exist close to 
lines with B(x) = Q @. 

For a pure el ectr ic modulation of the form V(x) — 
—Vocosqx, Eq. (Bl) is equivalent with a quadratic equa- 
tion for cosqx, which has a solution only if 1 — (qR) 2 + 
l(qR) 2 V /2E F } 2 > 0. This yields for the critical cyclotron 
radius 



2 = Vv_ 

qRcrit Ep 



l + y/l- (Vn/Ep) 2 



(B2) 



For Vb <C E F ,we recover the result B c 



B. 



Beton 



(V 



cqVo/ev F , which has previously been obtained by Be- 
ton et al. Q from the balance of the average Lorentz 
force and the maximum electrical force qV$ due to the 
modulatio n. F or stronger modulation, the B crit obtained 
from Eq. ( B2 ) increases faster than linearly with Vq and 
reaches the value B crlt = V2Bgf t on at V = E F . Then 
all trajectories are confined to a single period of the po- 
tential. For Bq > B cr n, i.e., for 2R < a\^2/n, all tra- 
jectories are "drifting" orbits, with left turning points on 
Xq{x) and right turning points on Xq(x), whereas for 
Bq < B cr it there exist also "channeled" trajectories hav- 
ing all their turning points either on x^(x) or on Xq(x). 
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For the step modulation V(x) = — Vq sign (cos qx), 
Xq (x) and Xq (x) are piecewise linear functions with dis- 
continuities at the walls of the potential wells, so that 
Eq. (Bl) can not be used to determine B cr n- Two types 
of channeled orbits may exist near a local minimum of 
Xq (x) (maximum of Xq(x)): traversing trajectories with 
turning points on opposite walls of a potential well, and 
skipping orbits which are reflected only by one of the 
walls. Traversing trajectories exist for sufficiently weak 
magnetic fields, B < B^j(Vb) = mcvp /eRf°£ with 



y/l + Vq/E F -y/1- V /E f 



(B3) 



Skipping orbits survive at arbitrarily high magnetic 
fields. Physically, this is a consequence of the fact that 
the singular electric field at a potential discontinuity 
never becomes negligible as compared to the Lorentz 
force due to Bq . For small modulation amplitude we find 
Bf r *%(V ) = (2/Tr)B Beton (V ). This indicates that the 
linear dependence of the critical magnetic field on the 
(weak) modulation strength holds independently of the 
detailed shape of the potential modulation, which may, 
however, affect the numerical prefactor. 

A characteristic feature of the channeled orbits is that, 
even at Bq ^ 0, their velocity in y-direction, v y = ui c (x — 
xq), does not change sign. 



2. Number of channeled orbits 

We now calculate the fraction of electrons occupying 
channeled orbits. The electron density n e i(x;E F ) — 
D [E F — V(x)]9(Ep — V(x)), i.e., the area density of 
occupied states at x with energy below Ep, is pro- 
portional to the area irv 2 (x) of the circle with radius 
v(x) = vp[l — V(x)/Ep] in v-space. The areal density 
of states at x with energy at E is given by the local 
density of states D(x;E) = dn e i(x; E) / dE — Dq. Both, 
n e i(x\ E) and D[x;E) are independent of the magnetic 
field Bq, and so are their spatial average values n e i{E) 
and D(E). An electron trajectory through x, at B = 0, 
is a channeled orbit, if the energy available for the mo- 
tion in the ^-direction is smaller than the height of the 
potential barrier, Ep — (m/2)vz < V max , where V max is 
the maximum value of V(x). With v y = v(x) sin ip, this 
means sin 2 <p > sin 2 ipo(x; Ep), where 



[x; Ep) = arcsin ^[E F - V max ]/[E F - V{x)} . (B4) 



Thus, the contribution of channeled orbits to the local 
density of states is D co (x;E) = Dq[1 — (2/n)(po{x;E)], 
and the density of electrons in channeled orbits (with 
energy below Ep) is n co (x;E F ) = D [E F - V(x)]0(E F — 
V(x)) if E F < V ma x, and 



,(x;E F ) 



2D r 



\jEp - V max \/V m ax ~ Vjx) 



(E F — V(x)) arcsin , 



'Vrr, 



V(x) 



E F - V(x) 



, (B5) 



if Ep > V max - For fixed average electron density n e i, the 
Fermi energy is independent of the modulation strength, 
if Vo < Ep = n e i/D . For V > Ep and B = 0, all tra- 
jectories become channeled orbits. For the step potential 
V(x) — — Vbsign(cosgx), channeled orbits exist only in 
the potential wells, and Eq. ( p5| ) yields the area density of 
channeled orbits in the well with V max — Vo and V[x) = 
— Vq. For Vo > Ep the requirement of fixed average den- 
sity fi e i yields a constant value of Ep + Vq = 2n e i/Do in 
the wells, so that for Vo > Ep the situation is identical 
to that at Vo = Ep, with vanishing electron density in 
the barrier regions. Since the local density of states is 
Dq, we normalize, for arbitrary Vo, the fraction of chan- 
neled orbits for the step modulation on the total density 
of states in the wells, ff o cp (B = 0, V ) = Df c P(x; E)/D 
(averaging also over the barrier regions, where no chan- 
neled orbits exist, would reduce f^ cp (B , Vq) by a factor 
of 2). For Bq = we obtain 



/ c s * ep (0,V ) = - arcsin J-^ 

7T V Ep 



Vo 



which for small e = Vo /Ep has the expansion 

^e P= 2^_(^)! + 0(e5/2)> 



(in 



(B6) 



(B7) 



This algebraic dependence on the modulation strength 
Vo is not an artifact of the step model. For the co- 
sine modulation V(x) = — Vocosgx one obtains, with 
f™(B ,V ) = D™(E F )/D = (D^(x,Ep))/D , 



r c ?(0,V ) = 2 - r-arcsinj^g-, (B8) 



E, 



V(x)> 



with the expansion 
4v/27 



£COS 

J CO 



7T" 



9tt 2 



0(e 



5/2\ 



(B9) 



for small amplitudes. 

For the step modulation, also the i?o-dependence of 
fco{Bo,Vo) can be calculated analytically. We suppress 
the lengthy result for < B < B%°£ [cf. Eq. @] and 
give only the result for B > Bf^, 



f 

■I ( 



step 



4i? 



V2e — \J\ — e arcsin 



2c 



(BIO) 



For the cosine model at < B < B™j t , we can cal- 
culate the positions and the values of the local extrema 
of Xq (x) [and similar for a;J(x)] analytically, and also 
D c c ° s (x; E). Contributions to the average of D c c ° s (x;E) 
come from the x-interval between the location x max 
of the local maximum and the value x upper satisfying 
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Xq upper) — %o {Xmax)- The value X U pp er and the inte 

gral of D c c ° s {x; E) between x max and x upper are calculated 
numerically. The results are plotted in Fig. ^| for both 
the cosine and the step model. Figure |^(a) shows the 
magnetic field dependence and (b) shows the dependence 
on the modulation strength of f C o(Bo,Vo), for Vq < Ep 
in all cases. For the cosine model at a given value of 
the magnetic field Bq, channeled orbits exist if Vo/Ep > 
(2/qR)[l - {qR)- 2 ] 1 / 2 . Their number increases with in- 
creasing Vq. For Vb — > Ep, fco(Bo,Vo) approaches, with 
vertical slope, a finite value, which decreases with in- 
creasing Bo. No channeled orbits survive if Bq becomes 
so large that 2/qR> \f2. For the step model, the given 
Bo-value sets a critical modulation strength V^ p given 
by V^f/Ep = (n/qR)[l - {n/2qR)*\V*. For V < V^g 
only skipping orbits exist, whereas for Vq > V^P also 
traversing orbits contribute to f C o(Bo,Vo), which leads 
to a change of the curvature of the curves in Fig. f7|(b) 
at Vq = V°r- t p . With increasing Vq, f C o(Bo,Vo) increases 
and reaches at Vq = Ep, also with vertical slope, the 
value 



f^(B Q ,E F ) 



arcsm 



b/y/2 



V2 i + yr 



b 2 /2 
(BH) 



if b = n/qR < V2, and f^(B ,E F ) = {2/n)V2/b if 
b > y/2. 

With the normalization chosen in Fig. 0(a), one ob- 
tains for the Bq dependence of f co nearly the same 
curve for all values of the modulation strength, < 
Vq < Ep. It can be shown analytically that, for 
weak modulation and to lowest order in Vo/Ep, a sin- 
gle curve should result for each of the models. Numer- 
ically we find this to be true for Vo/Ep < 0.2. For the 
step model / c s * cp « r - P (V ), V o )/f^(0, Vq) approaches 
the values 2/3 for Vq — ► and 2 /it for Vq — ► Ep. 

The algebraic dependence of the number of channeled 
orbits on the modulation strength leads to a similar al- 
gebraic Vq dependence of their contribution to the resis- 
tivity at Bq — 0, as we demonstrate now. 



3. Resistance at Bo = 

We now evaluate, for Bq = 0, the classical velocity- 
velocity-correlation formula for the conductivity in the 
extreme nonlocal limit Xq — ► oo. In this limit we expect 
the errors due to the violation of the continu ity equation 
to be small, as we discussed in section IV A 1 . In the limit 
-Bo = 0, the constant of motion is (v y ), the average veloc- 
ity in y-direction. The "drifting orbits" transform into 



orbits, which are unbound in x-direction and lead to a pe- 
riodic velocity v(x + a) = v(x), whereas the "channeled 
orbits" remain bound in ^-direction within one period, 
and their velocity in x-direction, v x , vanishes on average. 
Then the conductivity tensor reduces to |Q 



v?" = -inZ I —f ^v t ,((p,xo)v u (ip,x ) , (B12) 
v f/ 2 Jo a J~ir 2tt 

where V v {(p, Xq) is the average of a velocity component 
along the trajectory starting (for an arbitrary value of yo) 
at {xo,yo) with velocity v = vp{xo)(cos ip, sin cp), where 
vf(xq) =v F [l-V(x )/Ep} 1 / 2 . 

First, we consider a pure electric modulation, which 
implies that v y is constant along any trajectory. Since 
the modulation potential V(x) vanishes on average, we 



see immediately that a y 



do and <J X 



have the 



same Prud e valu es as the unmodulated system. Accord- 
ing to Eq. ( B12 ), the channeled orbits do not contribute 
to a xx . For the trajectories extended in x-direction, one 
has v x (ip, xq) = a/T(ip,xo), where 



T(<p,x ) 



dx 



o y /vl[l-V(x)/E F ]-v*(<p,x ) 



(B13) 



is the time an electron needs to traverse one period of 
the mo dulat ion. Since T{ip, xq) is an even f unctio n of 
ip, Eq. (B12) yields also a yx — 0. From Eqs. (B12) and 
(B13) we obtain for the electric cosine modulation 



-22- = (l + e)$, 



2c 



(B14) 



where e = Vo/Ep, 



I dt [t 2 K(t) VT^j}- 1 (B15) 

J n 



and K(t) is the complete elliptic integral of the first kind 
§§. Apparently a c ° x s -> for V — > Ep. For e -> 
one obtains a power expansion in yfk, the coefficients of 
which can be calculated numerically. To leading order in 
e one obtains 



3/2 



= 0.69454 (—) + O([V /Ep] 5/2 . 
Po \ E F/ 



(B16) 



For the step model the corresponding results can be 
calculated analytically. The exact result for the conduc- 
tivity is 



T stcp 

XX 

CTO 



1 



3- e 

7T V " 2e 
3 - 2e - 5e 2 



(B17) 



2c 



2ne ' V 1 + e 

with the leading term of the resistance correction 

— = ^{t L T + 0{[Vq/E f ]^). < B18 ) 
po 157r \Ep J 

Obviously, the step modulation yields the same half- 
integer power dependence of the resistivity on the modu- 
lation strength. Moreover, if we choose the mean values 
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of the squared modulation potentials equal, i.e., replace 
Vb i n the step potential by Vq/V2, the right hand side of 
Eq. QB18| ) becomes 0.571 (V Q /E F °f/ 2 where the prefactor 
is comparable with that in Eq. (|B16| ). 

The Bq — conductivities for pure magnetic modu- 
lations can be calculated similarly. For the modulation 
B rn (x) — B® n cos qx one obtains 



a 

w XX /-^ 

^0 



4o 



(B19) 



where $(77) is given by Eq. ( |B15[ ), and a = \/qR m = 
Wm/qvF- For a — > this yields a power expansion in 



similar to Eq. (B16) 



Po 



1.9645 (gi? m )- 3/2 + 0([gi? m ]- 2 ) , (B20) 



but now the correction term is larger. In this collisionless 
limit, the conductivity vanishes if all trajectories become 
bound in x-direction. For the cosine modulation this hap- 
pens at l/qR m = 1. For the magnetic step modulation, 
B m (x) — i?JJjSign(cos qx) the corresponding condition is 
a/AR m = 1, which means that a strip of constant B m (x) 
just can accommodate a cyclotron orbit with diameter 
2R m . As a function of a/AR m , the conductivity for the 
magnetic step modulation shows a similar behavior as 
that for the magnetic cosine modulation as a function of 
l/qR m - 

It can be shown analytically that, for Bq = 0, a yy — 
(To is not changed by the modulation. It is, however, 
interesting to calculate the fraction Act™ / cto contributed 
by the channeled orbits to a yy . This is easily shown to 
equal the fraction contributed to the density, 



Act co 

yy 

Co 



a dx n co (x; E F ) 
a fi e i 



(B21) 
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FIG. 1. (a) Sq ' ' vs. 1/qR oc Bo for different values of the 
mean free path A, given in terms of the modulation period 
a = 2ir/q; (b) and (c): selected magnetoresistance results 
for pure electric mod ulat ion Vix) — Vbcosga; to order V Q 2 , 
calculated from Eq. (4.1). In (b) Ap xx is plotted in units 
of p\ = (po/2)(Vo/Ef) , and the dash-dotted line represents 
the local limit for Xq — 1.0. In (c) Ap xx is plotted in units 
of p2 — Po{Vq I Ef) 2 (Xq /2j 2 , and the thin line represents the 
approximation of Ref. |LS|. 



FIG. 2. Magnetoresistivity vs. 1/qR oc Bo for (a) mixed 
electric and magnetic, and (b) mixed el ectr ic an d mobil- 
ity modulations, calculated from Eqs. (4.1)-(4.4) for co- 
sine modulations V(x) = Vbcosgx, Lo m (x) — lu^ cos qx, 
and r(x) = r^cosqx, of comparable effective amplitudes 
av = XqVo/2EF, as = tui^ and a M = rrj^, respectively, 
with Xq = 30. 7 is a small, but otherwise arbitrary, dimen- 
sionless constant. 



FIG. 3. Typical orbits (lower panel) in the electric su- 
perlattice potential V(x) — 0.4 Ef cos qx, with q — 2n/a, 
for three values of the magnetic field Bo corresponding to 
13 = 1/qR = 0.25, 0.1, and 0.001. The corresponding values 
of the constant of motion xo are indicated in the upper pan- 
els by horizontal bars ending on the curves x^ (x) defining the 
turning points. Channeled orbits exist for /3 > j3 cr it « 0.2. 



FIG. 4. Magnetoresistivity versus 1/qR oc Bo for a weak 
electric cosine modulation with Vo/Bf=0.02 and Xq=27.7. 
Solid lin e: n umerical result, dotted line: analytical approx- 
imation ( f4.l| ), dash-dotted: contribution of channeled orbits 
(see text). Inset: the same for magnetic cosine modulations 
with 1/qRm = 0.34, 0.15, and 0.05 (from top to bottom at 
arrow, for further explanations see text) and Ag=27.7. 



FIG. 5. Magnetoresistivity vs. 1/qR oc Bo for anisotropic 
scattering, modelled by Eq. (5.2), and an electric cosine mod- 
ulation with V /Ef = 0.02; r = 27.7 /qv F . 



FIG. 6. Magnetoresistivity vs. 1/qR for cosine mod- 
ulation of the electric potential and the mobility, with 
fvFq = 27.7. Solid lines: p xx /po for pure electric modulation 
with Vo/Ef ~ 0.2 and p yy /po for pure mobility modulation 
with r m = 1/r; broken lines: mixed modulation. 



FIG. 7. Fraction of channeled orbits, as function of mag- 
netic field B (a) and modulation strength Vo (b), for an elec- 
tric cosine modulation V(x) — — Vbcosga; and a step modu- 
lation V(x) — — Vo sign(cos qx); (a): B-dependence for fixed 
Vo/Ef = 0.01 (solid lines) and 0.95 (dash-dotted lines), (b): 
Vo-dependence for fixed 1/qR = 0.0, 0.05, 0.2, and 0.4 (from 
top to bottom), with solid lines for the cosine and dash-dotted 
lines for the step modulation. 
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